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Abstract. We identify regular structures in the globally chaotic spectra of an 
interacting bosonic quantum gas in tilted periodic potentials. The associated 
eigenstates exhibit strong localization properties on the lattice, and are dynamically 
robust against external perturbations. 
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1. Introduction 

Ultra-cold atoms in optical lattices provide a versatile testing ground for the physics of 
interacting many-body quantum systems, ranging from the characteristic properties 
of many-particle ground states P, over atomic transport [21 E], to the emergence 
of collective properties, and thermodynamic behavior [H |5]. The complexity of the 
many-body dynamics of such systems increases rapidly with the number of particles 
and lattice sites and has its characteristic spectral counterpart in a highly irregular 
parametric evolution of the energy levels |6l [71 [H |9] . This implies a high sensitivity of 
the time evolution with respect to changes in the initial conditions and/or perturbations 
of the generating Hamiltonian, and renders control of generic many-particle dynamics 
an extremely challenging task. However, nonlinear coupling can also give rise to the 
emergence of stable collective modes which opens new perspectives for robust control 
[T0[ [TT] . In our present contribution, we will identify such modes for the tilted Bose- 
Hubbard Hamiltonian (BHH), and demonstrate their pronounced localization properties 
as well as their extraordinary robustness against perturbations. In contrast to previous 
experimental and theoretical studies on bound-atom states (see, e.g., Refs. [121 [131 E]) 
the stability of these solutions is not a mere consequence of energy separation (by a 
spectral gap) due to the presence of interatomic interactions. Instead, we find bound 
states that exist within the bulk of the spectrum which can be considered chaotic in the 
sense of random matrix theory [T5] . 

2. Model 

The simplest quantum mechanical many-body description of ultra-cold bosons in a 
lattice is the Bose-Hubbard Hamiltonian [16] that incorporates both, the tunneling 
of individual particles between neighboring sites, and their pairwise on-site interaction. 
When the one- dimensional lattice is subject to an additional static tilt (due to, e.g., 
gravitation), the Hamiltonian takes the form 

J AI-l f/ ^ 

H = --Y^ {dl^di + h.c.) + - ^ n,(n, - 1) + F J] Ifii , (1) 
1=1 1=1 1=1 

where cii (a|) annihilates (creates) a particle in the Wannier state localized at the l-th 
site, fii = a\ai is the associated number operator, and M specifies the length of the 
lattice. Here we consider a tilt around the center of the lattice, and the on-site term 
/ hence takes the values /~ = -M/2 + / for even and [ = -(M + l)/2 + l for odd M. 
The BHH has two constants of motion, the energy E = (H) and the total number of 
particles N = The parameters J, U, and F describe the tunneling strength, 

the on-site interaction and the static tilting field, respectively. 

The model is based on a single (lowest) band approximation for the optical lattice 
[16] . This assumption is valid provided the kinetic energy, the interaction strength, 
and the local chemical potential (resulting from the tilt), are sufficiently small not to 
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Figure 1. (Color online) Spectrum of the Bose-Hubbard-Hamiltonian (Eq. ([T])) as a 
function of the static tilt F/J, with interaction strength U = I and tunneling coupling 
J = 1, for N = 3 particles and M — 11 lattice sites. The M eigenstates with the 
largest IPR (see Eq. Q) in the Pock number basis are plotted in red. a) A set of 
(almost) straight-line energy levels traverse the chaotic background, without changing 
the slope. This set is characterized by the localization of the corresponding solitonic 
eigenstates in the Fock basis, b) A zoom into the vicinity of the two lowest solitonic 
levels reveals small avoided crossings. 

excite higher Bloch bands. Therefore, the lattice needs to be sufficiently deep [T], 
to induce large band gaps, and the interaction energy must be smaller than the single- 
particle ground-state energy, in order not to modify the single-particle wave function 
considerably. In the experiment, these conditions can be met since all the parameters J, 
U, and F in the BHH are readily controlled While J and F are solely determined by 
the lattice geometry, the inter-atomic interaction U can additionally be adjusted using 
Feshbach resonances [TTl [18] . 

In the limit of large particle numbers ^ 1, the quantum dynamics may be 
described by the mean-field counterpart of the Bose-Hubbard Hamiltonian, the discrete 
Gross-Pitaevskii equation (GPE) (see, e.g., [191 l2D])- At fixed lattice length M, the 
mean-field limit is approached by increasing the particle number and, at the same time, 
keeping the scaled interaction UN at a constant value. As a mean-field approach, the 
discrete GPE does not explicitly contain the particle number N and, hence, hardly 
covers effects related to the granularity of matter [|} In the following, we are concerned 
with a rather small number of interacting bosons where such effects will be essential, 
and thus resort to the full many-body Hamiltonian ([T]). 

I We note that for small particle numbers, deviations of the GPE description from the many-body 
picture have been predicted, e.g., for a two-site lattice in Refs. [HIES, where a comparison to the 
Bose-Hubbard dynamics was performed. In Refs. [23l [21], both the BHH and the continuous GPE 
were further compared to the multiconfigurational time- dependent Hartree for bosons method, a more 
advanced propagation scheme for the many-body Schrodinger equation. The applicability of the GPE 
can though be extended, when combined with phase-space methods, i.e., when one considers not a single 
mean-field trajectory but propagates an ensemble of GPEs that reflects the initial quantum mechanical 
state (see, e.g., |22 and references therein). 
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Figure 2. (Color online) The average maximal population (AMP, red), Eq. ([s]), and 
the averaged inverse participation ratio (IPR, black), Eq. Q, of the eigenstates of the 
Bose-Hubbard Hamiltonian, Eq. ([T]). The solid lines correspond to the average of the 
M = 11 eigenstates with the largest individual IPR, the dashed lines to the average 
of all the other eigenstates. Within the range 0.1 ^ F/J < 0.4, the solitonic states 
are strongly localized, both on the lattice and in the Fock-space. The parameters are 
chosen as in Fig.jl] and the dimension of the Hilbert space is Af = 286. 




3. Parametric level evolution 

Tunneling, on the one hand, and interaction and tilt, on the other, define two 
(incompatible) symmetries of the system Hamiltonian ([T|. Hence, if either one of these 
terms dominates, the many-particle eigenstates exhibit the relatively simple structure 
of Bloch waves (for J > [/A^, F) or Wannier states (for UN > J, F or F > J, UN). In 
the generic case, however, when all three terms have comparable weight, good quantum 
numbers with an unambiguous labeling of the system eigenstates cannot be defined, 
since the energy levels exhibit a complicated parametric evolution with U, J, or F, and 
avoided crossings of variable size abound |8l [9l EU ET] . This is the spectral manifestation 
of quantum chaos, and is nicely illustrated by the overwhelming number of energy levels 
in Fig.[TJ 

However, the figure also shows energy levels with constant slope over a wide interval 
of the static tilt F (at fixed U and J), which thus represent many-particle eigenstates 
with characteristic properties which are invariant under changes of F: By virtue of the 
Hellmann-Feynman theorem [28] , 

where (•) represents the expectation value with respect to a single such state, the slope of 
these energy levels defines a constant center of mass of the many-particle wave function, 
and, hence, suggests invariant localization properties of the particles on the lattice, under 
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changes of F. This is corroborated by the average maximal particle number (AMP) of 
these states on a single lattice site 

AMP = ^^^^^^^iM^^ , (3) 

in Fig.[2| in comparison to the AMP of the chaotic, i.e., irregular states of the spectrum. 
The average (■) in ^ runs over the respective sample of states, at fixed value of F. 
Clearly, in the irregular eigenstates the particles are largely spread out over the lattice, 
with AMP < 0.25, while the solutions with constant slope in Fig.[l]are stronger localized 
on the lattice. Whereas this localization is to some extent present already for vanishing 
tilt, it becomes strongest once the corresponding energy levels cross the bulk, i.e., for tilt 
strengths 0.1 < F/J < 0.4, where more than 60% of all {N = 3) particles are localized 
on a single site. Because of their pronounced localization properties, that we will show 
to remain unchanged under variation of the tilt, we refer to these as solitonic states in 
the following [29] . In contrast to previously identified localized states [IS, E] , the 
presently discussed solitonic states are not merely a consequence of energetic isolation 
with respect to the remaining part of the spectrum. Finally, we note that the observed 
localization is consistent with the fact that there exist as many solitonic states as sites 
of the lattice - M = 11 in the case investigated in Figs. [l]and[2] 

Given the fact that the solitonic states are localized around one lattice site, we 
further check whether they also exhibit localization in Fock space. A measure for 
the latter is the inverse participation ratio (IPR) with respect to the Fock basis 
\bj) = |nj^i,nj 2, ■■■Uj^M)-, which forms the eigenbasis of the BHH for vanishing tunneling 
J = 0. The IPR is defined as: 

M 

IPR(|^)) = 5^|c,r , (4) 

i=i 

where the Cj are the expansion coefficients of the state in a given basis, and M is the 
Hilbert-space dimension. The IPR represents the inverse number of basis states that 
are occupied by the state and varies from unity - when \%p) coincides with one basis 
state - to l/M, for lip) an equally-weighted superposition of all basis states. In fact, 
the states highlighted in red in Fig. [T] are precisely those eleven states with the largest 
IPR, i.e., those which exhibit the strongest localization properties in Fock space. To 
quantify this statement we calculate the IPR averaged over the eleven solitonic states, 
as well as for all non-solitonic states [§J For the latter. Fig. [2] shows a moderate and 
monotonic increase of the IPR with growing F, what can be attributed to the onset 
of Stark localization, as expected for a very large tilt (see, e.g. In contrast, the 

IPR of the solitonic states is about one order of magnitude larger than the IPR of the 
irregular states, for 0.15 ^ F/J < 0.35. For stronger tilt, some of the sohtonic states 
disappear what we discuss in Section Generating mechanism further down. 

§ The fluctuations of AMP and averaged IPR for tlie solitonic states, as visible in Fig. [2] can be 
attributed to avoided crossings that locally affect the levels and are not washed out due to the rather 
small number of solitonic levels that enter the average. For the non-solitonic states, AMP and averaged 
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Figure 3. (Color online) Two spectra of the Bose-Hubbard-Hamiltonian (Eq. ([T])) as 
a function of the static tilt F/J (cf. Fig.[l]). The system parameters are a) interaction 
strength U — 0.5, tunneling coupling J = 1, TV = 6 particles, M — 7 lattice sites, with 
Hilbert-space dimension Af — 1716, and h) U — 0.3, J — 1, N ^ 10, M ~ 5, and 
Af ~ 1001. In both panels, the M eigenstates with the largest IPR (see Eq. Q) in the 
Fock number basis are plotted in red. The insets show the spectrum on a larger scale. 

A few words regarding the dependence of our observations on the system parameters 
are in order. As we will explain in detail in Section 5, the presence of solitonic states 
is not linked to a specific set of parameter values, but we expect them quite generally 
to appear in a regime where the Bose-Hubbard spectrum can be considered primarily 
chaotic. This expectation is corroborated by Fig. |3} which shows the parametric 
evolution of two Bose-Hubbard lattices with a) = 6 particles in M = 7 lattices 
sites, and b) A^ = 10 particles in M = 5 lattice sites, and interaction values U chosen 
in the chaotic regime [27j. In agreement with Fig. [l| we observe as many solitonic 
states (marked in red) as sites of the lattice - M = 7 and M = 5 in Fig. [sja) and b), 
respectively. In the following, wefixf/ = J = l,M = ll sites, and A^ = 3 bosons, as a 
typical case of well-developed quantum chaos in coexistence with solitonic states. 



IPR show less fluctuations, since the average is taken over significantly more states. 
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Figure 4. (Color online) Left: Magnification of the spectrum where the energy levels 
for the initial states |^) are marked. Black circles/red lines depict irregular/solitonic 
initial states. Right: Dynamical evolution of the average IPR, Eq. Q, in the 
instantaneous basis for various rates R, Eq. ([s]), with Fi/J = 0.1 and Ff/J = 0.4. 
Solid/dashed lines represent solitonic/irregular initial states. 



4. Dynamical stability 

The solitonic states' invariance properties spelled out by the constant slope of their 
energy levels have an additional expression in the comparatively small avoided crossings 
with the irregular states of the spectrum, as evident from Fig.[T]D). This implies small 
coupling matrix elements between solitonic and irregular eigenstates, and therefore 
suggests an enhanced stability of the solitonic states under variations of F. 

To test this conjecture, we investigate the multi-particle dynamics under a linear 
ramp of the tilt, F{t) = Fi + Rt, from the initial value F{0) = Fi, to the final value 
F{At) = Ff. That is, the tilt is varied with the slew rate 
Ff - F 

for both solitonic and irregular eigenstates |\E') of the BHH chosen from the same energy 
range in the bulk of the spectrum, as shown in Fig.|4^). The stability of the initially 
prepared states is characterized in terms of the time-dependent IPR (see Eq. (|4])), for a 
given rate R. 

Let us first consider the IPR in the instantaneous basis, i.e., in the Hamiltonian's 
time- dependent eigenbasis that follows the evolution of the static field F{t). It is 
shown in Fig.|4]D) and reveals the broadening of the initial state due to transitions 
to other modes. We observe a strikingly different behavior for both initial conditions: 
Whereas the IPR of the irregular states decays rapidly, the solitonic states' IPR decreases 
significantly slower. We also find that the decay rate observed in Fig. [Ijo), is essentially 
independent of R, after an initial transient. A complementary view on the stability can 
be obtained from the IPR in the initial (fixed) basis defined by Fi, see Fig. [5^). For very 
large rates, R = 1, the change in the tilt F occurs on a time scale much faster than the 
internal (tunneling) dynamics on the lattice. Hence, in this diabatic regime the system 
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Figure 5. (Color online) Left: Dynamical evolution of the average IPR, Eq. in the 
initial (fixed) cigenbasis for various rates R, Eq. with Fi/J = 0.1 and Ff/J = 0.4. 
Solid/dashed lines represent solitonic/irregular initial states. Right: Drop of the IPR 
against the rescaled time i = t^fR. For short times, the curves for all three rates fall 
on top of each other, thereby forming the upper group of irregular states and the lower 
group of solitonic states. That is, both types of initial conditions exhibit excellent 
agreement with linear response theory - which predicts |31j a super-Gaussian decay 
for the survival probability P{t) ^ exp[— (the dash-dotted line has slope four, and 
is drawn to guide the eye). Note, however, the offset between irregular and solitonic 
states which amounts to about one order of magnitude. 



wave function cannot adapt to the changing potentiaL Consequently, the IPR remains 
largely unaffected and a difference between solitonic and irregular states is hardly visible. 
In contrast, as we reduce the ramp, we observe more and more pronounced drops of the 
IPR for both types of initial states. However, the drop in the IPR of the irregular states 
is about one order of magnitude larger than that of the solitonic ones, what spells out 
the stability of the latter. 

We finally remark that the IPR's sensitive dependence on the slew rate i?, present 
in the instantaneous basis shown in Fig.|4|3), can be understood from linear response 
theory. To this end, we note that, for sufficiently short times, the inverse participation 
ratio IPR(|\E'(t))) = \ (\Efj| tj (t) l^^o) I'' is essentially given by the square of the survival 
probability P{t) = \ (^o| U{t) |^o) P since | U{t) |^o) P ^ for j ^ 0, where V{t) 
is the time-evolution operator. For linearly driven chaotic systems, quantum linear 
response theory predicts a super-Gaussian decay P(t) ~ exp {-RH'^] [211, what suggests 
a scale invariance of Pit) and hence of the IPR with respect to the scaled time t = ty/R. 
This is confirmed in Fig.[5|D) where for solitonic and irregular initial states the same 
dependence is observed. We stress, however, that between the two groups of curves 
there is an offset of about one order of magnitude, that is, the absolute decay is much 
weaker for the solitonic states ^ 

II Note that the latter exhibit oscillations of 1 — IPR around t = 1, which are likely to result from the 
regularity of the integrable part of the spectrum. 
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5. Generating mechanism 

In the previous section, we numerically confirmed the dynamical stability of the solitonic 
states, as suggested by their parametric level evolution. We now turn to the discussion 
of the underlying mechanism |32] ■ The existence of eigenstates of a many-body system 
with all particles localized close to each other, despite the presence of repulsive inter- 
particle interactions U, was experimentally first demonstrated for pairs of atoms [12]. 
Theoretical investigations of bound atom states as well as of bound quasi-particle states 
in nonlinear lattice systems, (see, e.g., [33l |3ll |35l |36l |371 |38l |39]) also revealed the 
existence of bound states including more than two (quasi-) particles [131 [HI HOI SH 112] ■ 
Recently, also Bloch oscillations for initially localized states of interacting bosons 
were investigated [IS]- Such repulsively-bound many-particle states are formed as a 
consequence of the energy mismatch between the on-site interaction energy and the 
maximal kinetic energy that can be realized in the lowest energy band in the lattice. 
That is, they are energetically isolated from the remainder of the spectrum. 

In contrast to that, for our presently solitonic states, which run through the bulk of 
the spectrum, also energetically allowed transitions (to the irregular states) are blocked, 
and, in addition, lead to an enhanced stability under perturbations. This can be 
understood in terms of the transitions between the energy eigenstates that are induced 
by the tilt, as discussed in the following. The time-independent part of the corresponding 
transition amplitude from some initial state to a final state |\E'/) is given by the 
corresponding matrix element of the center-of-mass operator (\E'j-| Ifii l^'j). As shown 
by our analysis in Fig.[2] above, in case of an irregular state the atoms are distributed 
over the entire lattice while the solitonic states are distinguished by the fact that most of 
the atoms occupy the same lattice site. As a result, the solitonic states are approximate 
eigenstates of the center-of-mass operator and hence the matrix element (\E'/| Ifii |\E'j) 
with f) an irregular and |\I'j) a solitonic state becomes very small, since both are 
system eigenstates and thus, mutually orthogonal. More intuitively, a transition from a 
solitonic to an irregular state requires the redistribution of essentially all atoms over 
the entire lattice. The tunneling of an atom over more than a single lattice site, 
as well as the simultaneous tunneling of more than one atom, is, however, negligible 
for typical lattice-depths, since it corresponds to higher-order processes in the Bose- 
Hubbard Hamiltonian. Whereas such processes can occur as resonance phenomena for 
large interaction strengths U [13] , they are negligible in the presently discussed regime of 
intermediate interactions. Furthermore, in a tilted, infinitely long lattice (but with finite 
particle numbers), all states are localized [13] and in that case, even analytic estimates of 
the localization volume of the states (and thus of the above transition matrix elements) 
can be obtained via composite-particle eigenvectors (see [13]). 

One might now wonder whether, similarly to the solitonic states, there are system 
eigenstates where (to a good approximation) all atoms but one are located on a single 
site, and the remaining particle is localized on some other site. Indeed, we find such 
states, that show comparable features as the solitonic states described so far, and they 



Robust states of ultra-cold bosons in tilted optical lattices 



10 



will be referred to as solitonic states of second order. For the case of three atoms, such 
states have also been analyzed in Ref . [H] . For a transition from a solitonic state of first 
order to a solitonic state of second order only tunneling of a single atom over a single 
site is necessary, and hence one expects significantly larger coupling as compared to 
irregular states. Thus, in the simplest approximation one can treat two such states as 
a two-level system of (idealized) Fock states |. . . , 0, A^, 0, . . .) and |. . . , 1, — 1, 0, . . .). 
Assuming vanishing coupling J = 0, their energy would coincide only at F = U{N — 1) 
which, for [/ = 1, is far outside the regime in which solitonic states exist (see Figsjl] 
and|2|. Once we take into account the tunneling, an avoided crossing of width 2J\fN 
emerges between the two states and thus affects the stability of the solitonic states for 
tilts beyond the threshold value Ft ^U{N — V) — J\fN . For the parameter values used 
in Figsjl] and |2| this corresponds to Ft ~ 0.27, and this is indeed where the solitonic 
states start to dissolve. 

If, for larger iV, more particles are bound in a solitonic state, more have to undergo 
tunneling processes in order to transform into an irregular state. From this point of view, 
the stability of the solitonic states is expected to be enhanced with increasing particle 
number at constant interaction strength U . Yet, for an increasing number of particles 
located on a single site, three-body interactions that are not accounted for by the BHH 
might become non-negligible. These three-body collisions may result in the formation 
of untrapped molecules leading to additional decay channels [1] for solitonic states. It is 
therefore crucial that the condition UN = const, is met by the experiment, since then, 
three-body processes are not enhanced due to an effective reduction of the interaction 
strength U d& N is increased. The same condition is essential for the construction of 
a mean-field phase space for the Bose-Hubbard Hamiltonian, for large boson numbers 
and constant lattice length M (see, e.g., P, [191 HH HHl HE] ) : The underlying principle is 
that under the requirement UN = const, the mean-field dynamics is not changed as 
is increased, due to the scaling behavior of the BHH The phase space will, in general, 
be mixed regular and chaotic [19]. The solitonic states should then be identifiable with 
regular islands that correspond to highly localized mean-field solutions as the ones that 
have been found, e.g., for three site systems (see, e.g.. Fig. 6 of |16] and inset of Fig. 6 
in |17]). Semiclassical arguments then guarantee the existence of quantum mechanical 
eigenstates which are localized within these islands and are only weakly coupled to the 
states living on the chaotic sea. 

While the analysis of the mean-field limit is beyond the scope of the present paper, 
we would like to shortly address some related aspects. Due to the comparably small 
particle number < 10 in the spectra presented in Figs. [T] and [3} these systems are still 
well below the mean-field limit. Instead, we show in Fig. [6] how the parametric level 

% We note that the parametric level dynamics presented in Figs, [l] and [s] obey this requirement, i.e., 
UN = 3 in the three cases shown. The mean-field dynamics is not the same, however, among lattices of 
different size M: Given the fact that the lattice lengths M were on purpose chosen differently, such as 
to demonstrate the general existence of solitonic states for various configurations, the respective phase 
spaces have different dimension of 2M [12] and are thus not identical. 
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dynamics itself evoWes as the particle number is increased, for fixed lattice size M and 
constant effective interaction UN = 3. Owing to numerical manageability, we chose a 
three-site lattice (M = 3), the smallest possible system for which chaos can be expected 
[19j . and increase the particle number from = 10 to iV = 50. We first remark that 
in all cases solitonic states are indeed observed: We highlighted the six levels with the 
highest IPR and note that in all three panels, the red level that crosses the bulk with 
negative slope and that is marked in red up to around F/J^ 0.6, corresponds to one 
of the three solitonic states of first order. Second, qualitatively speaking, larger particle 
numbers result in a smaller effective Planck constant: That is, for increasing N and 
UN = const., more quantum states share the available, iV-independent phase space. 
Thus, also more states can exist within the island, i.e., an increased number of solitonic 
states of higher order is expected. This, as well, is confirmed by Fig. |6| For = 10 
particles, we find essentially one solitonic state with negative slope, corresponding to 
localization on the first lattice site, with negative single-particle on-site energy Fl = —F, 
see below Eq. ([T]). In contrast, for = 20, two such levels can be clearly distinguished 
while for A^ = 50 already three such states are among the six levels with the highest IPR. 
Finally, the detailed inspection of the spectra shown in Fig. |6] reveals that the critical 
tilt strength Ft, at which the solitonic state of highest order disappears, displays a weak 
dependence on the particle number A^. Namely, for A^ = 10 (A^ = 50) we determine 
Ft to be approximately F/J ^ 0.65 {F/J ^ 0.75). This indicates that for increasing 
A^, the criterion for Ft should be based on the following semiclassical consideration 
instead of the simplified picture given above: As the tilt F is increased, the stable 
island (associated with localization on the corresponding lattice site) shrinks and, thus, 
fewer and fewer solitonic states can reside on it. Finally, at a certain value of the tilt 
F, given by the (A^-independent) mean-field equations, the island vanishes and this 
determines the critical value Ft at which the solitonic state of highest order disappears 


With respect to the experimental viability of the investigated setup, it is fair to say 
that not only the lattice geometry is under exquisite control by the experimentalists, 
but also interaction strength and particle detection: In typical optical lattice-based 
experiments, the ratio U/J can be sensitively adjusted by a variation of the lattice 
depth [IHl [5]. Depending on the atomic species, U can additionally be controlled, 
independently of the coupling tunneling J, via Feshbach resonances [171 [IB]- For 
example, in experiments with Cs atoms, the scattering length could be tuned from 
small negative to very large positive values, over several orders of magnitude [51]. In 
the same setup, magnetic levitation was rapidly switched off within about hundred /xs 
and was used to compensate for gravitation. Presumably, this levitation can as well be 
tuned to a nonzero value as to adjust the strength of the gravitational forces, what would 
amount to controlling the lattice tilt F. Finally, we note that the detection efficiency has 

+ For a two-site Bose-Hubbard system a similar strategy was successfully applied in Ref. [35] (see also 
references therein), and a tilted chaotic three- mode system has been studied in the mean- field limit in 
Refs. [HEn]. 
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Figure 6. (Color online) The spectrum of the Bose-Hubbard Hamiltonian (Eq. ([!])) 
for fixed effective interaction UN = 3 and three different particle numbers N, plotted 
as a function of the static tilt F/J (cf. Fig.[T]). For better comparability the energies are 
additionally scaled by N. The system parameters are tunneling coupling J = 1, Af = 3 
lattice sites, and particle numbers (Hilbert-space dimension) a) = 10 {Af = 66), b) 
N = 20 {Af = 231), and c) N = 50 {Af = 1326). In aU panels, only the six eigenstates 
with the largest IPR (see Eq. Q) in the Fock number basis are plotted in red to keep 
the graphical presentation clear. We note, though, that the number of solitonic states 
increases with N. That level which crosses the bulk with negative slope and, in all 
three panels, is marked in red up to around F/ J k, 0.6 corresponds to a solitonic state 
of first order. 



reached the single-atom level [521 [53], what permits a very precise counting of occupation 
numbers. Thus, the experimental detection of solitonic states should be within reach. 

6. Conclusion 

In summary, we have studied the dynamical evolution of ultracold bosons confined on a 
one-dimensional optical lattice, which is subject to a tilt. Our present work was focused 
on the parameter regime in which kinetic, inter-particle interaction, and on-site energy 
are balanced. In this regime, where the corresponding Bose-Hubbard Hamiltonian with 
an additional tilt term can be characterized as being primarily chaotic in a spectral 
sense [SU [261 [27], we have identified regular structures in the parametric level evolution. 
Associated with these structures are the solitonic states, which are distinguished by a 
strong localization on the lattice as well as in the Fock space, which barely changes with 
the tilt strength. Unlike other studies on energetically isolated bound states jUl |13l [H] , 
this property results from the solitonic states' weak coupling to the bulk of energy levels 
[32]. 

Based on the inverse participation ratio, evaluated in the fixed as well as in 
the adiabatic basis of the Bose-Hubbard Hamiltonian, we compared the dynamical 
stability of these solitonic states to an ensemble of neighboring states in the bulk of the 
spectrum, as the tilt was ramped with different slew rates. We observed a drastically 
increased robustness of the solitonic states, spelled out by a significantly larger inverse 
participation ratio at the end of the dynamical evolution. 

The present work belongs to that branch of ultracold atom physics, which is 
concerned with comparatively small particle numbers. This direction has recently seen 



Robust states of ultra-cold bosons in tilted optical lattices 



13 



substantial experimental advances, not at last in the resolution of imaging techniques 
[52t |53] that would allow to directly identify the solitonic states. Their remarkable 
dynamical robustness makes them excellent candidates, e.g., for the preparation of stable 
quantum few-body states, in a parameter regime where, due to the generic presence of 
spectral chaos, a substantial portion of the eigenstates sensitively depends on system 
parameters such as the tilt. On the other hand, provided that coherent superpositions 
of solitonic states can be prepared, reasonably stable tunnel dynamics of more than one 
boson [13] would be realizable for relatively small interaction values. Finally, in the 
light of the mean-field limit discussed in the previous section, the crossover in stability 
of these states from the few- to the many-body regime might be explored. 
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